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Abstract. We give the map representing the evolution of a qubit under the action of 
non-dissipative random external fields. From this map we construct the corresponding 
master equation that in turn allows us to phcnomenologically introduce population 
damping of the qubit system. We then compare, in this system, the time-regions 
when non-Markovianity is present on the basis of different criteria both for the non- 
dissipative and dissipative case. We show that the adopted criteria agree both in the 
non-dissipative case and in the presence of population damping. 



Comparison of non-Markovianity criteria in a qubit under random external fields 2 

1. Introduction 

In quantum systems the dynamics of decoherence, and that of quantum correlations, 
is qualitatively different if the environment is Markovian (without memory) or non- 
Markovian (with memory) [TJ [2j |3]. For example, for composite quantum systems 
independent non-Markovian environments, entanglement may present revivals [U El E] 
or trapping [7J Ej defending it against sudden death [9]. Non-Markovian systems are 
utilized in several physical contexts such as quantum optics [?], solid-state physics 
|10j . quantum chemistry [11] and quantum information processing [12] . It is therefore 
essential to establish criteria to identify and quantify the non-Markovian behavior in an 
open quantum system. Among the criteria, one introduced by Breuer-Laine-Piilo (BLP), 
is based on the concept of temporary flow of information from the environment back into 
the system and quantify non-Markovianity as an increase in the distinguishability of two 
evolving quantum states [13]. A second one, due to Rivas-Huelga-Plenio (RHP), instead 
measures the deviation of the dynamical map from divisibility [H]. A third one has 
been also proposed by Andersson-Cresser-Hall (ACH) that uses the negative decoherence 
rates, appearing in the master equation, as a primary measure to completely characterize 
non-Markovianity [15] . An all-optical experiment has been recently developed to control 
transitions from Markovian to non-Markovian dynamics [16] . 

A natural question is then if the different criteria agree in identifying non-Markovian 
behaviors in the system dynamics. It has been shown that, for a qubit coupled 
to environments via the Jaynes-Cummings or dephasing models, the BLP and RHP 
criteria have exactly the same non-Markovian time-evolution intervals and therefore are 
equivalent [17] . In an analysis performed for a driven qubit in a structured environment 
it has been suggested that the two measures may disagree [18] and successively it has 
been shown both for a classical and a quantum toy-model [3]. Comparisons among 
the three criteria, including the ACH one, showing possible non-equivalence in realistic 
systems are instead still missing. 

In this paper we address this issue. In particular our aim is to verify, for a realistic 
physical system made of a qubit subject to random external fields both with and without 
dissipation, if the BLP, RHP and ACH criteria give concordant results in individuating 
the non-Markovian time-regions in the system dynamics. 

2. Model 

We consider a realistic system made of a qubit subject to random external fields both in 
a non-dissipative and in a dissipative case. In the following we describe the two cases. 

2.1. Non-dissipative random external fields 

Our system is a qubit interacting with an environment constituted by a classical field 
mode with fixed amplitude but with random phase equal either to zero or to it with 
probability p — 1/2. This model has been introduced to study the possibility of revivals 
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of quantum correlations in absence of back-action [T5] and describes a special case of 
a qubit subject to a phase noisy laser [201 EI]- The dynamical map is of the random 
external fields type [22| 123] and, in the qubit basis {|1), |2)}, is written as [19] 

i 2 

A(t,0)p(0) = -J2u i (t)p(WHt), (1) 

i=i 

where 

JT( .,( cos(At) e-^sm(Xt)\ 
i[) [-e^sm(Xt) cos(At) J' [) 

with i — 1,2 and 0i = 0, 02 = 7r. = e~ tHit / h is the time-evolution operator 

associated to the Hamiltonian = ih\(a + e~' l ^ i — a-e 1 ^*), where cr+, cr_ are the qubit 
raising and lowering operators and A is the qubit-field coupling constant that depends 
on the field amplitude. The Hamiltonian Hi is given in the interaction picture (rotating 
frame) at the qubit-field resonant frequency u. 

In order to use the non-Markovianity measures introduced above, the knowledge of 
both dynamical map and master equation is required. In our model we directly have 
the map and we also have to construct the corresponding master equation. To obtain 
the master equation starting from the map of Eq. (CD) we follow the procedure proposed 
in Ref. [20] which gives (details of calculations are reported in the Appendix) 

dp/dr = Lp{r) = tan2r(crj,pcr y — p), (3) 

where r = Xt is a dimensionless time. It is worth to observe that this form of master 
equation, associated to our system, presents a time-dependent rate, tan(2r), which is 
the same that has been previously introduced only formally in a general master equation 
to study non-Markovian behavior [T3l [H] . 

2.2. Dissipative case 

The model of random external fields described above is non-dissipative and can be 
generalized to a dissipative case. Although it is not obvious to introduce a source 
of dissipation directly into the map, it is simple to do into the master equation. We 
phenomenologically add population damping with rate 7, in the standard Lindblad form 
with generator 70". [23] , into the master equation of Eq. ([3]) which now becomes 

dp/dr = Lp{r) = tan2T(a y pa y — p) 

+ 7(a_pa+ - pa + (j_/2 - a + a_p/2), (4) 

where 7 = 7/ A is a dimensionless decay rate. In the following, we shall use the map of 
Eq. ([T]) and the master equations of Eqs. ([3]) and (j4]) to analyze if the different criteria 
individuate the same time-regions when non-Markovian behavior occurs. 

3. Comparison among the criteria in the non-dissipative case 

We shall first apply the three non-Markovianity criteria (BLP, RHP and ACH) to the 
case of non-dissipative random external fields. 
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3.1. BLP criterion 

The BLP criterion is based on the distinguishability of two evolving quantum states 
quantified by the trace distance [13], that is D(pi(t), P2(t)) = \\\pi(t) — p 2 (t)\\i where 
= TrV A^A, pi{t) = A(t,0)pi (i = 1,2), whose variation rate is 

a(t)=dD( Pl (t),p 2 (t))/dt. (5) 

The dynamical map A(t, 0) is non-Markovian, according to BLP, if there exists a pair 
of initial states p±, p 2 such that for some time t > the distinguishability of the two 
states increases, that is a(t) > 0. This is interpreted as a flow of information from the 
environment back to the system which enhances the possibility of distinguishing the two 
states. 

Let us apply this criterion to the model of non-dissipative random external fields. 
Chosen two arbitrary initial states 

( u) ae^\ ( p (3e^\ 

Pl= {ae^ 1- U )> P2= [pe-^ 1 - M J ' (6) 

and substituting them into Eq. (jHJ) we obtain 

tr(r) = -VSsin(4r)/|6|, (7) 

where a = (p — u) 2 + (a cos <pi — j3 cos ip 2 ) 2 and b = cos 2 2r + a sin ipi — f3 sin ip 2 - The 
sign of this quantity does not depend on the value of the parameters of the initial 
states and thus permits a general comparison with the other criteria. In particular it 
is readily found that <t(t) > (i.e., the dynamics exhibits non-Markovianity) when 
7r/4 + k(ir/2) < r < (k + l)vr/2, where k is a non-negative integer number. 

3.2. RHP criterion 

The RHP criterion is based on the divisibility of a dynamical map and it is independent 
of the system state. If the map A(t, 0) is divisible, it satisfies the condition A(t+ e ,o) = 
A( t+e)t )A( t) o) (e is a time interval) that is usually attributed to Markovian evolution. It 
is possible to show that the map A(t, 0) is completely positive (CP), and then divisible, 
if and only if (A t+e t ® 1 2 ) | *^ > ) | > 0, where |$) is a maximally entangled state of 
two qubits (one is subject to the map while the other is the isolated ancilla) and 11 2 
is the two-dimensional identity matrix [T4]. For a qubit subject to a master equation 
dp/dt = L t (p), where L t is a Lindblad operator, in the limit of e — >• the solution 
(dynamical map) of this equation formally tends to A t + e ,t ~ > eLte - Expanding this 
solution up to the first order in e it is possible to introduce the quantity [1] 



(() = Um II [l 4 + e (L8 1 ! ,)]|»)<»||| 1 -l 

e^0+ e 

where ||v4||i indicates the trace norm. It is shown that g(t) > if and only if the original 
map A(t, 0) is indivisible, that is exhibits non-Markovian behavior. 

In our case of non-dissipative random external fields, identifying L t with that of 
the master equation of Eq. ([2]), we obtain g(r) = — 2tan2r if tan2r < and g{r) = 
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otherwise. It is immediately seen that non-Mar kovian behavior (g(r) > 0) occurs 
just in the same temporal regions individuated above by the BLP criterion, that is 
vr/4 + k(n/2) < r < (k + 1)tt/2. 



3.3. ACH criterion 

This criterion is based on the property of complete positivity (divisibility) of the 
dynamical map deduced through the sign of time- dependent decoherence rates that 
may appear in the master equation. This criterion is also independent of the system 
state. Consider a qubit governed by a master equation in the canonical (Lindblad-type) 
form, in the interaction picture [15] 

% = E^^(t)p4(t) - l -L\{r)L k {r)p - l -pL\{r)L k {r)l (9) 

k 

where the traceless operators L k (r), time-dependent in general, describe different 
decoherence channels and 7fc(r) are the corresponding decay rates that can be also time- 
dependent. The different decay channels are orthogonal in the sense that Tr(LtL fc ) = 5 jk . 
If the 7fc(r) are positive at all times, then the time evolution is completely positive in 
any time interval with a Mar kovian behavior. On the other hand, if some of the ~f k (t) is 
negative, the time evolution exhibits non-Markovian behavior that can be then naturally 
characterized by the function f k (r) = min^/^r), 0] for each decoherence channel |15j . 
This criterion is conceptually similar to the RHP one and it is convenient due to its 
immediate application once having the expression of the master equation. 

In the master equation of Eq. ([3]), associated to our model of a qubit under non- 
dissipative random external fields, the only (dimensionless) decay rate is tan2r. Once 
again we find that the time regions when non-Markovian behavior occurs correspond to 
the negative values of tan 2r. 

The above results show agreement among the three criteria in individuating time- 
regions of non-Markovianity here considered, in the case of non-dissipative random 
external fields. 



4. Comparison among the criteria in the dissipative case 



We now analyze RHP and ACH criteria in the case of a qubit subject to random external 
fields and to population decay, whose master equation is given in Eq. (J4]). We do 
not consider the BLP criterion that requires the knowledge of the qubit evolution and 
therefore the solutions of the master equation of Eq. (jlj): this will be treated elsewhere. 
The function g(t) of Eq. ([8]) of the RHP criterion now becomes 

-7/2 - 7 i(r)/2 + (V2/4) [g + {r) + g_{r)\ , (10) 



9\T) 



where g±(r) = < 7 2 + [7 + 71 (r)] 71 (r) ± a/7 2 + 7i( r ) f an d 7i( r ) = 2tan2r 



1/2 
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Figure 1. Comparison between the function g(r) of the RHP criterion (blue solid 
line) and the function /(r) of the ACH criterion (red dashed line) as a function of the 
dimensionless time r, for a dimensionless decay rate 7 = 3. There is non-Markovianity 
when g(r) > according to RHP and when /(t) < according to ACH. 



To use the ACH criterium, we put the master equation of Eq. (@| into the canonical 
form of Eq. ([9]) by using the procedure of Ref . [15] ; two orthogonal decay channels arise 
with rates 

7 ± (r) = (7 + 2 tan 2r ± V7 2 + 4 tan 2 2r)/2, (11) 

and corresponding operators L± = X)i=i 2 ^i a i/V^, where o"j (i = 1,2) are the usual 
Pauli matrices and 

( ±) _ i(-2 tan 2r ± ^ + 4 tan 2 2r) 

y7 2 + (2 tan 2r =f ^T 2 + 4 tan 2 2r) 2 

f/ 2 (±) = 7/y / 7 2 + (2tan2r^ v / ^ 2 + 4tan 2 2r) 2 (12) 

Being 7-(r) < 7+(t) at any time, the non-Markovianity regions according to ACH are 
characterized only by the function /_(r) = min[7__(r), 0]. From Eq. (1111) . the condition 
7-(t) < is satisfied when 47tan2r < (i.e., 7r/4 + k{n/2) < r < (k + l)n/2). 
Therefore, the ACH criterion in the dissipative case individuates non-Markovianity in 
the same time regions of the previous non-dissipative case. 

In this dissipative case, the ACH criterion evidence non-Markovian behavior in the 
same time regions individuated by the RHP criterion. This is displayed in Fig. [H where 
it is seen that the function g{r) of the RHP criterion is greater than zero exactly when 
the function /(r) of the ACH criterion is lower than zero. All the above results are 
independent of the initial state of the system. 

5. Conclusions 

In this paper we have analysed three different criteria (BLP, RHP and ACH) identifying 
non-Markovian behaviors in a realistic system made of a qubit subject to random 
external fields, both in a non-dissipative and in a dissipative evolution. We have first 
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exactly obtained the master equation corresponding to the qubit dynamical map of 
random external fields. We point out that the form of the master equation, associated to 
our system, contains the time-dependent rate tan(2r) that has been previously inserted 
only formally in a general master equation to study non-Markovian behavior [131 E3] • 
We have then phenomenologically introduced population damping directly in the master 
equation associated to the map of random external fields. 

We have found, in the non-dissipative case, that the three criteria agree into 
individuating non-Markovianity time-regions. For the model of random external fields 
with population decay, both the RHP and ACH criteria individuate the same time 
regions of non-Markovian behavior. 

The results of this paper may provide new insight to the topic of characterizing the 
non-Markovianity in a realistic open quantum system. 
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Appendix A. Master equation associated to the model of non-dissipative 
random external fields 

In this Appendix we summarize the steps to obtain the master equation of Eq. ([3]) 
from the map of random external fields of Eq. (CQ) by following the general procedure 
described in Ref. [20] . 

The general steps are as follows. Let us apply a map A(t, 0) to the basis operators 
Gi = Oij 'a/2 (i = 0, . . . , 3), where a® = 1 and the remaining o~i are the Pauli matrices, 
and define a matrix F with elements = Tr[GfcA(t, 0)(Gi)]. The idea is of constructing 
a matrix FF^ 1 (or, more generally, FF if F is not invertible). In our case F is invertible 
and it is possible to calculate the matrix R, with elements defined by 

Rab = ^(FF-^MGrT^n], (A.l) 
rs 

where r a = |ai)(a2|, U = \Pi){P2\, with |ai), |a 2 ) and being the qubit basis 

states |1), 1 2). The general expression of the master equation is then 

L(p(r)) = p{r) = Rabit)r a p{t)rl (A.2) 

ab 

where the operators r are To = |2)(2| = o" + cr_, t\ = |1)(1| = o"~cr+, r 2 = |2)(1| = a + and 
r 3 = |1)(2| = cr_, with cr± = (01 ± i£7 2 )/2. In our case of random external fields with 
the map given in Eq. (pQ), we obtain the matrix F 

( 1 \ 

cos2r 



10 
cos2r j 



(A.3) 
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from which one easily obtains the matrices F 1 , F and therefore the matrix FF 1 . 
Chosen the basis { 1 2) (2 1, |1)(1|, |2)(1|, |1)(2|} and using Eq. flXj} , we find the R matrix 



as 



R 



I 



\ 



tan 2r 
tan 2r 







tan2r 
tan2r 










tan2r 
- tan 2r 






- tan 2r 
tan 2r 



/ 



Finally, using Eq. f ]A.2j) we obtain the desired master equation 
dp/dr = tan 2r{a y pa y — p). 



(A.4) 



(A.5) 
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